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ABSTRACT 


In  this  report  we  present^  method  for  calculating  the  aerodynamic 
influence  coefficients  (AICs)  based  on  third-order  piston  theory  with  an 
optional  correction  to  agree  with  Van  Dyke's  quasi-steady  second-order 
theory.  The  AICs  are  computed  assuming  the  airfoil  to  have  a  rigid  chord 
with  or  without  a  (rigid  chord)  control  surface.  The  influence  coefficients 
relate  the  surface  deflections  to  the  aerodynamic  forcesU^rough  the  following 
definitions .  Tn  the  oscillatory  case, 


=  pa)Vs[C^]  Jj^- 


and  in  the  steady  case, 


/ 


fJ  =  (l/2)pv2(S/c)  [Cj^J  {h| 


The  piston  theory  is  limited  to  high  Mach  number  (or  high  reduced  frequency) 
but  Van  Dyke's  quasi-steady  correction  extends  the  validity  to  some  lower 
supersonic  Mach  number  at  low  reduced  frequency, 

/^he  Aerospace  IBM  7090  Computer  Program  Number  HMll  provides 
the  AICs  from  this  theory  in  both  a  printed  and  an  optional  punched-card 
output  format.  The  program  capacity  is  25  surface  strips,  15  Mach  numbers 
and  20  reduced  velocities  for  each  Mach  number./^  V  ^ 
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SECTION  I 


FORMULATION  OF  THE  PROBLEM 


A.  Introduction 

The  pressure  on  a  lifting  surface  is  normally  given  by  a  surface 
functional  relationship.  However,  in  the  limit  of  a  high  Mach  number  (or 
high  reduced  frequency),  this  relationship  becomes  a  point  function.  As  a 
consequence  of  this  limit,  aerodynamic  influence  coefficients  (AICs)  may  be 
speciTied  exactly  by  a  strip  theory,  and  control  surface  and  camber  effects 
may  be  determined  in  a  straightforward  manner. 

The  present  formulation  derives  the  AICs  from  third-order  piston 
theory  for  a  lifting  surface  with  control  surface  (both  assumed  rigid  in  the 
chordwise  direction;  i.  e.  ,  no  camber  is  presently  considered).  The  deriva¬ 
tion  differs  only  slightly  from  that  of  Ashley  and  Zartarian^  in  that  in  the 
present  case  the  third-order  pressure  coefficient  is  generalized  to  account 

for  sweep  and  steady  angle  of  attack,  and,  following  a  suggestion  of  Morgan, 

2 

Huckel,  and  Runyan,  a  correction  (optional)  is  suggested  to  give  agreement 

3 

with  the  second-order  quasi-steady  supersonic  theory  of  Van  Dyke.  This 
quasi- steady  correction  should  extend  the  validity  of  the  piston  theory  to 
lower  supersonic  Mach  numbers  at  low  reduced  frequencies.  The  derivation 
given  here  is  taken  from  Ref.  4;  further,  the  computational  aspects  of  the 
present  report  are  an  extension  of  the  computing  procedure  of  Ref.  4. 

B.  Sign  Convention 

The  flutter  sign  convention  is  used  in  the  oscillatory  case:  forces  and 
deflections  are  positive  down;  rotations  are  positive  with  leading  edge  up. 

The  aerodynamic  sign  convention  is  used  in  the  steady  case:  forces  and 
deflections  are  positive  up;  rotations  are  positive  with  leading  edge  up. 
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C.  Derivation  of  Equations 


We  quote  here  the  development  of  Miles  in  obtaining  the  piston  theory 
pressure  coefficient.  There  are  two  cases  of  interest.  Th^ firsj/ assumes 
that  the  angle  of  attack  is  small  enough  that  there  are  pressu^ifperturbations 
on  the  expansion  side  of  the  surface.  The^^^s^cor^assumes  that  the  angle  of 
attack  is  large,  and  that  the  expansion  pressure  approaches  a  vacuum  and  is 
ineffective  in  producing  perturbations.  Because  of  the  difficulty  in  specifying 
the  transition  from  low  to  high  angle  of  attack,  we  shall  restrict  the  present 
consideration  to  the  first  case,  the  low  angle  of  attack. 


"Hayes'  hypersonic  approximation  states  that  any  plane  slab  of  fluid 
initially  perpendicular  to  the  undisturbed  flow  may  be  assumed  to  remain  so  J 
as  it  is  swept  downstream  and  to  move  in  its  own  plane  under  the  laws  of 
one -dimensional,  unsteady  motion.  Thus,  the  probleni  of  a  wing  having  an 
a rbitrarily  prescribed  mntinr)  nornral  to  its  surface  may  be  reduced  to  the 
considerati on  of  the  _p ne jr-dimfijggipnal  motion  of  a  piston  into  an  otherwis e 
undisturbed  flaw.  This  problem  is  relatively  simple  if  the  disturbances 
produced  by  the  piston  are  treated  as  simple  waves,  for  then  the  pressure 
on  the  piston  depends  only  on  the  instantaneous  velocity  there,  w\.  and  is 
given  by  » 

p/p  =  [  1  +  (1/2)  (y-l)  (w/a  )] 


tj/ 

''fv.o 


(1) 


where  p^  and  a^  are  the  values  of  pressure  and  sonic  velocity  in  the 
undisturbed  flow. 

"The  result,  Eq.  (1),  is  exact  for  an  expansion,  but  the  presence  of  a 
shock  front  (and  consequent  departure  from  isentropic  flow)  renders  it  only  ap¬ 
proximate  for  a  compression.  Lighthill  has  suggested  a  cubic  approximation 
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to  be  adequate  fox*  pi*actical  application  if  |w/a^|  <  1.  The  aeries  expansion 
yields 


p/p^  =  1  +  y(w/a^)  +  (1/4)  y(v  +  1)  (w/a^)^ 
+  (1/  12)  v(y  +  1)  (w/a^)^ 


(2) 


Lighthill  has  shown  that  this  expression,  Eq.  (2),  is  within  six  percent  of  the 

value  given  by  either  Eq.  (1)  or  the  exact  solution  with  the  shock  at  maximum 

5 

permissible  strength.  " 


The  pressure  coefficient 

2 

noting  that  q  =  {.y/2)  p  M  . 


C 
.  P 


=  (p  -  Pq)/*!  is  found  from  Eq.  (2)  after 


Cp  =  (2/M^)  [(w/a^)  +  (1/4)  (y  +  1)  (w/a^)^ 

-1  (1/12)  (V  -t  1)  (w/a^)^]  (3) 

Following  a  suggestion  of  Morgan,  Huckel,  and  Runyan,  ^  we  may  generalize 
this  result,  Eq.  (3),  by  writing 

Cp  =  (2/M^)  (Cj(w/a^)  +  C2(w/a^)^  +  C3(w/a^)^l  (4) 

in  which  for  piston  theory 

Cj  =  1,  C2  =  (y+1)/4.  C3  =  (y  +  1)/12  (5) 

and  for  the  quasi-steady  theory  of  Van  Dyke^ 

Cj  =  M/p,  =  [m‘*(Y,+  l)-4p^] /4p'*,  C3  =  (y+1)/12  (6) 


3 


Van  Dyke  gives  only  the  second-order  solution  so  that  the  value  of  is 
taken  from  the  piston  theory  result.  The  use  of  the  modified  coefficients 
and  could  extend  the  low.er  Mach  number  limit  of  piston  theory. 

We  may  now  calculate  the  lifting  pressure  coefficient  from  Eq.  (4)  and 
the  local  piston  velocity.  The  normal  velocity  (positive  away  from  the 
surface)  on  the  upper  and  lower  surfaces  of  a  symmetrical  thin  airfoil  having 
thickness  distribution  2g(x)  and  angle  of  attack  is  given  by 

w  =  V  (g  -  Q  -  v)  ,  (7a) 

u  '®x  o' 

Wjj  =  V  (g^  +  +  v)  ,  (7b) 

where  v  is  the  unsteady  component  of  the  dimensionless  downwash. 

For  the  case  of  small  angles  of  attack,  the  lifting  pressure  (positive 
down)  is 

=  =  -  (4/M)  [(C  +  2C  Mg 

P  P^  Pi  1  Z  X 

+  3C,M^g^)  (q  +v)  +  C-M^(a  +v)^]  .  (8) 

3  ®x  o  3  o  ■' 

If,  consistent  with  the  small  perturbation  assumptions  of  aeroelastic  analysis, 
only  the  terms  linear  in  v  are  retained,  Eq.  (8)  becomes 

Cp  =  -  (4v/M)  [Cj  +  ZC^Mg^  +  3C3M^(g^  +  q^)]  (9) 

Before  discussing  the  swept  wing  transformation,  it  is  appropriate  to 

review  the  limitations  of  Eq.  (9).  Ashley  and  Zartarian^  have  shown  that  the 

2 

piston  theory  is  applicable  if  any  of  the  conditions  M  >  >  1,  Mk  >  >  1,  or 


4 


k  >  >  1  is  met.  We  see  that  for  low  reduced  frequency  the  Mach  number 
necessarily  must  be  high.  However,  if  the  reduced  frequency  is  large  the 
Mach  number  is  not  necessarily  large;  in  fact  it  could  be  transonic  or  even 
subsonic.  At  this  point  it  is  apparent  that  any  sweep  correction  introduced 
to  bring  piston  theory  into  line  with  linearized  supersonic  theory  must  be 
considered  as  a  low  frequency  approximation. 

The  result,  Eq.  (9),  applies  to  the  swept  wing  case  if  all  quantities  are 
determined  by  the  flow  characteristics  normal  to  the  leading  edge.  The 
expressions  may  be  rewritten  in  the  form 

(Tp  =  .  (4v/M)  +  3C3M^gi  +  Q^)]  (10) 

The  transformation  from  the  normal  values  to  the  free  stream  values  are  the 
following : 


the  Mach  number 


the  geometry 


=  M  cos  A 

X  =  X  cos  A 
K  =  b  cos  A 


the  angles  of  attack  and  slope 


a  =  a  /cos  A 
o  o 


P  =  p/  cos  A 
®x  ^  S^/cos  A  ; 


the  dynamic  pressure 


q  =  q  cos  A 


(lla) 

(llb) 

(llc) 


(lid) 

(lie) 
(Ilf) 


dig) 
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and  the  pressui’e  coefficient 


C  =  C  cos^  A  .  (1 ih) 

P  P 

We  note  that  h  and  k  are  invariant.  From  the  dimensionless  downwash 

V  =  (l/V)|h  +  Va  +  (x  - 

+  [Vp  +  (x  -  Xj^)p]  Ux  -  ^ 

whicl)  for  harmonic  motion  becomes 

V  =  ikh/b  +  [  1  +  i(k/b)  (x  -  Q 

+  {  1  +  i(k/b)  (x  -  Xj^)]  P  J.(x  -  Xj^)  ,  (13) 

we  find  the  transformed  value 

V  =  ikh/^  +  [  1  +  i(k/S)  (x  -  x^)]'ia 

+  [  1  +  i(k/b)  (x  -  x^)  ]  p  _l(x  -  Xj^)  ( 14a) 

=  ikh/b  cos  A  +  [  1  +  i(k/b)  (x  -  x^)]  o-/ cos  A 

+  [  1  +  i(k/b)  (x  -  Xj^)]  (p/cos  A)  ^^(x  -  Xj^)  {14b) 

=  v/ cos  A  (14c) 
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The  transformed  pressure  coefficient  becomes 

2  2 

C  =  "C  cos  A  =  [  -4(v/cos  A)  cos  A  /M  cos  A  1 
P  P  . 

X  [  C,  +  2C_(M  cos  A  )  (g  / cos  A  ) 

lb  X 

+  3^2{M  cos  a  (g^  +  Q^)/ cos^  A  ]  ( 15a) 

=  -(4v/M)  +  Z'C^Mg^  +  3C2M^(g^  +  a^)  ]  .  (15b) 

We  note  that  the  sweep  effect  shows  up  only  in  the  coefficients  and  ^2' 
for  piston  theory,  there  is  no  effect 

^^=Cj  =  l,  ?T2  =  C2  =  (y  +  l)/4  ,  (16) 

and  for  the  quasi-steady  supersonic  theory 

TTj  =  M/(M^  -  sec^A  )^^^ 

^2  =  [M^y  +  U-4  8ec^A(M^  -  sec^  A)] /[4(M^  -  sec^  A)^  (17) 

Equation  (17)  is  seen  to  be  the  most  general  result.  If  sec  A  is  taken 
as  zero  then  the  piston  theory  results,  Eqs.  (16),  are  obtained;  and  if  sec  A 
is  taken  as  unity  the  sweep  correction  is  not  made  in  the  quasi-steady  super¬ 
sonic  result. 

We  next  consider  the  integration  of  the  pressure  coefficients  obtained 
above.  The  oscillatory  aerodynamic  coefficients  referred  to  the  leading  edge 
are  defined  by  the  following  equations. 
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(18a) 


dL/dy  =  4p<jJ^  ^  ^  ^3  ^ 

o  o 


dM/dy  =  4pw^  b^/M^  h^/b  +  ^3  ^ 

O  0^0 


(18b) 


dT/dy  =  4pu^  b^(T^  h^/b  +  T^  .  a  +  T  p 

O  0  *^0 


(18c) 


The  lift,  moment,  and  hinge  moment  are  found  from  the  pressure  coefficient 


2b 


f  CO 

dL/dy  =  q  J  dx 


(19a) 


2b 


^  cu 

dM/  dy  =  q  J  x  dx 


(19b) 


r2b 

dT/dy  =  q  j  (x  -  Xj^)  C  dx 


h'  p 


(19c) 


where  the  pressure  coefficient  is  given  by  Eq.  (15b). 


Cp  =  -(4/M)  [):,  +  4  3C,  M^g2  4  »^)  I 

X  jikh^/b  +  [  1  +  ikx/b]  a  +  [l  +  ik(x  -  Xj^)/b]p_l(x  -  Xj^)  j  (20) 


and  we  have  taken  the  pitch  axis  at  the  leading  edge  x^  =  0.  We  define  the 
following  dimensionless  thickness  integrals 


rZh 

Ij  =  (l/2b)  j  g^dx 


(21a) 
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I 


=  (l/4b^) 

rZh 

1  *  Sx 

(21b) 

=  (l/8b^) 

f^'’  2  . 

X  g  dx 

•^0 

(21c) 

=  (l/2b)  ^ 

0 

(21d) 

=  (l/4b^) 

z  , 

1  X  g  dx 

•^0 

(21e) 

=  (l/8b^) 

Z  Z  - 

(21f) 

=  (1/Zb) 

-Zb 

/  Kx'** 

(ZZa) 

-Zb 

J,  =  (l/4b^)  f  X  g 


dx 


(ZZb) 


Zb 


■»  r  z 

Jj  =  (l/8b  )  J  X  8x 


(ZZc) 


Zb 


=  (1/Zb)  J  dx 


(ZZd) 
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(22e) 


,2b 


=  (l/4b  )  j  X  dx 


=  (l/8b^)  f  dx  .  (22f) 

These  thickness  integrals  are  evaluated  at  the  end  of  this  section  for  a 
typical  airfoil.  If  we  substitute  Eq.  (20)  into  Eqs.  (19)i  make  use  of  the 
definitions  Eqs.  (21)  and  (22)  of  the  thickness  integrals,  and  identify  the 
resulting  expressions  with  Eqs.  (18),  we  obtain  the  oscillatory  aerodynamic 
coefficients 


=  -  iK  /k 

(23a) 

o 

L 

n 

=  -  K  /k^  -  iK  /k 

(23b) 

o 

=  -  K^/k^  -  MKg  -  2K4|^)/k 

(23c) 

=  -  iK  /k 

(23d) 

o 

M 

n 

=  -  K  /k^  -  iK-/k 

(23e) 

o 

“Po 

=  -  -  i(K^  -  2K5ei^)/k 

(23f) 

=  -  i  (Kg  -  2K^  e^)/k 

(23g) 

T 

a 

o 

=  -^^5  ■  -  i(K^  -  2K5ej^)/k 

(23h) 
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I 


T  =  -  (Kg  -  2K^eh)/k^  -  i(K^  -  +  4K^d)/k  (23i) 

where 

4j^  =  Xh/2b  (24a) 

Kj  =  (1/M)  ["CT^  +  Z'C^MIj  +  +  Q^)  ]  (24b) 

=  (1/M)  +  4'C^Ul^  +  3CgM^(2l5  +  a^)]  (24c) 

Kg  =  (4/3M)  +  bTT^MIg  +  3CgM^(3I^  +  q^)]  (24d) 

=  (1/M)  ^j(l  -  ej,)  +  2^2^^!  ■*■  ^CgM^L 
Kg  =  (1/M)  Cj(l  -  ej)  +  4Z^IA3^  +  3CgM^[2Jg  +  a^(  1  -  i^)]  (24f) 

Kg  =  (4/3M)  T:^(1  -  e‘^)  +  6C2MJ3  +  3C3M^[3Jg  +  a^(l  -  e^)]  .  (24g) 


To  conclude  the  derivation  of  the  oscillatory  aerodynamic  coefficients, 

we  calculate  the  thickness  integrals  for  the  typical  airfoil  of  Fig.  1.  We 

approximate  the  airfoil  by  two  parabolas  and  a  line.  The  equation  of  the 

•S' 

forward  parabola  that  goes  through  the  leading  edge  and  is  horizontal  at 
the  point  of  the  maximum  thickness  is 

gj(x)/c  =  (t/2)  (x/x^)  (2  -  x/x^)  (25) 

The  approximation  by  a  sharp  leading  edge  is  consistent  with  the  theory 
having  ruled  out  detached  shock  waves. 
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where  t  =  t  /c.  The  second  parabola,  horizontal  at  the  point  of  maximum 
max 

thickness  and  going  through  the  hinge  line,  is 


g^(x)fc  =  (t/2)  jl  -  (1  -  rj^)  [{x  -  (26) 

where  r,  =  Tv/t  and  t,  =  t, /c.  The  line  connecting  the  hinge  line  and  blunt 
h  n.  n  h. 

trailing  edge  is  given  by 


g^(x)/c  =  (tj^/2)  [  1  -  (1  -  r^)  (x  -  x^)/(c  -  x^)]  (21) 
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where  ~  differentiating  we  find  the  desired  slopes 


g^(x)/c  =  (t/x^)  (1  -  x/x^) 

(28a) 

g'{x)/c  =  -  t(1  -  r,  )  (x  -  X  )/(x,  -  x 
**2  h  m  '  h  m 

(28b) 

g^(x)/c  =  -  {tj^/2)  (1  -  r^)/(c  '  Xj^) 

(28c) 

From  the  slopes,  the  thickness  integrals  follow  immediately, 
control  surface  integrals  first  yields 

Computing  the 

^1  = 

^2=  . 

f  g^de  =  -  (1/2)  (tj^  -  T^) 

^h 

/‘eg5de  =  -(i/4)(r^- Vd  +  ep 

^h 

(29a) 

(29b) 

^3'-^ 

j'  -(i/6)(T^-T,)(i  +  e^  +  e^) 

^h 

(29c) 

II 

I  =  (i/4)(tj^  -  T^)^/{i  -  ej,) 

^h 

(29d) 

f  eg|de  =  ( 1/8)  ( 1  +  ^)/(  1  -  i^) 

^h 

(29e) 

'  r  Inside  =  (I/IZ)  (T^  -  t/  (1  +  +  Ih)/(1  -  e^) 

^h 


(29f) 


The  complete  airfoil  integrals  Vjecome 


I 


I 


I 


r  ^ 

2  "  I  +  ^2  =  -  ■*■  +  ^2 

2 

3  =  i  e'g^de  +  ,t3 

=  (t/12)6^  -  (1/12)  (r  -  T^)  (31^  +  26^5^  3  i^)  +  J, 

4  =  jT "  g|<ie  +  J4  =  .  ( 1/3)  (T  -  T^)"/(e^  -  e„)  +  j 
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(30a) 


(30b) 


{30c) 


(30d) 


=  t^/12  +  (1/12)  (t  -  Tj^)^  (3|j^  +  en,)/(eh  *  +  J5  (30e) 

i6  =  {  ”  +  j, 

=  (-"/30)e^  +  (1/30)  (T  -  T^)^  (65^  +  3e^5„,  +  -  5^)  +  j^  (300 

where  ^  =  x/c,  £  =  x  /c,  and  £,  =  x. /c. 

TTt  m  n  n 
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Having  obtained  the  oscillatory  aerodynamic  coefficients,  we  are 
now  in  a  position  to  derive  the  AICs.  We  consider  the  given  and  equivalent 
force  systems  in  Fig.  2.  The  equivalent  forces  are  arbitrarily  placed  at 
the  quarter-chord,  the  control  surface  hinge  line,  and  the  trailing  edge. 
The  derivation  must  relate  the  forces  F^,  F^,  F^  to  the  deflections  h^, 
h.^^,  h^  through  the  given  leading  edge  aerodynamic  coefficients  and  deflec¬ 
tions  h^,  a,  (3.  We  begin  with  the  force  equivalence. 


1 

1 

1  ^ 

■^r 

L 

o 

b/2 

{b/2  +  d) 

2b 

^2 

,  =  < 

M 

o 

0 

0 

c 

a 

T 

o' 

(M) 


The  loads  and  deflections  are  related  through  the  definitions  of  the  oscillatory 
coefficients. 


r  L  'I 

1  0  0 

o 

M 

o 

►  =  4pu)^b^Ay 

0  b  0 

T 

0  0b 

1  o  J 

The  equivalence  in  the  deflections  is  given  by 


o 

(1  +  b/  2d) 

-b/2d 

0 

’^l' 

ba 

.  = 

-b/ d 

b/d 

0 

^2 

be 

V  J 

b/d 

-{b/d  +  b/ c  ) 
a 

b/  c 

.^3  . 

(33) 
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Fig.  2.  Original  (a)  and  Equivalent  (b)  Force 
Systems  and  Geometry  for  Oscillatory  Case. 
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Substituting  Eq.  (33)  into  (3Z),  Eq.  (32)  into  (31),  and  solving  for  the  forces  yields 


(1  +  b/2d) 

-b/d 

)  (3b/2d  -  1) 

^2 

►  =  4pc 

0  b  Ay 

-b/2d 

b/d  -{b/ 

c^) (3b/2d) 

^3 

0  • 

0 

b/  c 

a 

— 

L  L„ 

n  R 

r 

(1  +  b/2d) 

-b/2d 

0 

o 

O 

X 

M  M. 

-b/d 

b/d 

0 

o 

T  T. 

a.  Lj 

b/d 

-(b/ d  +  b/c  )  b/ c 

o 

o 

J 

From  the  definition  of  the  AIC  matrix 


f|  =  p(uVs  [C^]  |h 


(34) 


(35) 


and  by  identity  with  Eq.  (34),  we  find  the  AICs  for  a  single  strip. 


(1  +  b/2d) 

-b/d 

(b/c^) (3b/2d  -  1) 

4(b/b^)^(Ay/s) 

-b/2d 

b/d 

-(b/c J  (3b/2d) 
a 

0 

0 

b/ c 

a 

o 

0 

d 

(1  +  b/2d) 

-b/2d 

0 

X 

u 

n 

“b 

-b/d 

b/d 

0 

o 

vX 

o 

o 

T 

a 

o 

b/d 

-(b/d  +  b/c^) 

b/  c 

a 

(36) 
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In  the  absence  of  a  control  surface  Eq.{36)  reduces  to 


[C^]  =  4(b/b^)^  (Ay/s) 


'(1  +  b/2d) 
.  -b/2d 


(37) 


The  complete  AIC  matrix  for  a  surface  of  N  strips  appears  in  the  partitioned 
form 


in  which  the  first  null  partition  is  reserved  for  control  points  at  which  the 
aerodynamic  forces  are  negligible  (e.  g.  ,  external  stores)  and  in  which  the 
remaining  partitions  are  of  the  size  3  x  3  or  2x2  according  to  whether  or 
not  the  strip  has  a  control  surface. 
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The  steady  AIC  matrix  follows  from  the  oscillatory  solution  as  a 
limiting  case.  If  we  compare  the  definition  of  the  steady  matrix 

|f|  =  (l/2)p  V^S/c)  [C^J  |h|  (39) 

with  the  oscillatory  definition  Eq.  (35),  we  observe 

[  C,  ]  =  2(sc/S)  lim  k^[  C,  ]  (40) 

hs  r  h 

r 

From  the  previous  section  we  find  the  limiting  values  of  the  oscillatory 
coefficients  to  be 


lim 

k  -0  \  '’o 

r 


lim 

k^-*0 


k^L 

r  a 


o 


lim  k^M 
1  rv  r  a 
k  -*0  o 


r 


lim  k^T 

1  n 

k  -*0  o 


lim 

k  -0 
r 


,  k^M,  .  k^T,  \  =  0 
r  h  ’  r  h  / 
o  o/ 

=  -  Kj(b^/b)^ 

=  -  K^(h^/h)^ 

=  -  (K5  -  2K4ei,)  (b^/b)^ 
=  -  K4(b^/b)^ 


(41) 


(42a) 


(42b) 


(42c) 


(43a) 
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(43b) 


lim  k^M  =  -  K  (b  /b)^ 
k  -*0  ^  Pq  ^ 

r 


lim  k^T  =  -  (Kg  -  (b^/b)^  (43c) 

k  -►0 
r 
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SECTION  II 


GENERAL  DESCRIPTION  OF  INPUT 


A.  Units 

Since  all  dimensional  input  is  geometrical  and  the  aerodynamic  matrix 
is  dimensionless,  only  a  consistent  set  of  length  units  is  neces sary--inches  or 
feet. 


B.  Classes  of  Numerical  Data  and  Limitations 

The  data  required  by  the  program  are  control  and  option  indicators, 
geometry,  Mach  numbers,  and  a  set  of  reduced  velocities  for  each  Mach 
number.  The  example  problem  illustrates  their  use. 

1 .  Example  Problem 

We  consider  the  four-strip  wing  shown  in  Fig.  3  at  Mach  numbers 
1.  8  and  2.  5.  We  use  reduced  velocities  of  4.  0  and  8.  0  for  both  Mach  num¬ 
bers,  and  compute  the  steady  case  for  Mach  2.  5.  The  aerodynamic  matrices 
will  be  computed  by  piston  theory  and  by  Van  Dyke's  quasi- steady  variation. 
Strips  2  and  3  are  considered  to  have  control  surfaces.  The  thickness 
integrals  will  be  computed  for  an  assumed  airfoil  (constant  across  the  span) 
having  10  percent  thickness,  maximum  thickness  at  40  percent  chord,  and  a 
blunt  trailing  edge  having  1.  5  percent  thickness. 


2.  Program  Restrictions  and  Options 


a.  The  number  of  strips  into  which  a  wing  may  be  subdivided 


must  be  <  25. 

b. 


The  number  of  Mach  numbers  must  be  <15. 


c.  The  number  of  reduced  velocities  used  with  any  one  Mach 
number  must  be  <  20. 
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Strip  No. 

Ay(ft) 

b(ft) 

Calft) 

d(ft) 

1 

4.  7 

12.  28120 

0 

11.  9 

2 

4.  2 

9.  50000 

5.  25000 

9.  0 

3 

3.  6 

7.  06250 

3.  99375 

6.  6 

4 

3.  1 

4.  96875 

0 

4.  5 

Strip  No. 

m 

^h 

T 

^h 

* 

1 

0.  4 

(not  used) 

0. 

1 

0.  015 

(not  used) 

2 

0.  4 

0. 72368421 

0. 

1 

0.  050 

0.  015 

3 

0.  4 

0. 71725664 

0. 

1 

0.  050 

0.  015 

4 

0.  4 

(not  used) 

0. 

1 

0.  015 

(not  used) 

sec  A  =  1 . 

25 

s  = 

554.  0 

sq  ft 

b  =6. 
r 

5  ft 

c  = 

21.  0 

ft 

s  =  15 

.  6  ft 

a  's 
o 

(cons 

tant)  =  5. 

0° 

N.B.  The  trailing  edge  thickness,  is  listed  as  the  hinge  line  thickness 
in  the  case  of  no  control  surface. 
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d.  If  it  is  desired  to  compute  the  steady  matrix  [  ^  zero 

or  negative  value  of  V/b^co  must  be  supplied  to  the  program.  (S  and  c  must 
also  be  provided.  ) 

e.  Thickness  integrals  may  be  given  or  computed.  If  given 

they  may  be  given  only  once  with  each  deck  and  are  considered  constant  with 
strips,  may  be  constant  or  vary  with  strips  (for  each  Mach  number). 

(t,  t^)'s  maybe  constant  or  vary  with  strips.  and  maybe 

constant  or  vary  with  strips. 

f.  The  control  surface  strips  must  be  a  continuation  of  the 
main  surface  strips;  e.  g.  ,  in  the  case  of  a  partial  span  control  surface  the 
inboard  and  outboard  span  stations  should  be  used  as  boundaries  of  the  main 
surface  strips. 

g.  As  many  complete  sets  (decks)  of  input  data  may  be 
supplied  as  desired  (one  following  the  other). 


23 


SECTION  III 


DATA  DECK  SETUP 


A.  Loading  Order 

Input  decks  punched  from  keypunch  forms  are  loaded  behind  column 
binary  deck  HMll.  The  data  for  each  deck  should  be  in  the  following  order: 

(1)  Heading  Card  1 

(2)  Heading  Card  2 

* 

(3)  NTHRY,  NTHICK,  NALPHA,  NTAUS,  NZETAS 

(4)  ISZ,  MSZ,  NOPUNJ,  JSZj,  JSZ., . -^^^IvISZ 

(5)  sec  A,  b^,  s,  S,  c 

(6)  Ayj.  Ay^ . Ajg^ 

^2 . ^ISZ 

^al’  ^a2 . ^alSZ 

^2 . '^IGZ 

(10)  Machj,  Mach2,  ....  Machj^g^; 

(lla)  If  thickness  integrals  are  given: 

(a)  When  all  c^^  =  0  tabulate  only  Ij,  l^,  I^. 

(b)  Any  c^j.  f  0  then  include  Jj,  J^i  •  .  .  i  •  i  and 

^hl'  C'  •  •  •  •  ^hKZ  »'  =  1  only  i.  needed). 

(llb)  If  thickness  integrals  are  computed: 


(a)  T,,  Tj^j,  r^.  ;  Tjg^,  Tj^jg^,  r^jg^ 

[  if  NTAUS  =  1  only  Tj,  and  are  needed;  if  c^^  =  0  (i.  e.  , 
no  control  surface),  the  trailing  edge  thickness  (t^.  )  is  listed  as 
Tj^j,  and  the  location  for  may  be  left  blank  for  these  strips]. 


Please,  no  remarks  about  our  Greek! 
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(12a) 


^ml’  ^hl’  ^m2’  ^h2’  •  '  ’  '  ^mlSZ  '  ^hlSZ  NZETA 

only  i  ,  and  ,  are  needed;  if  c  .  =  0,  the  program 
^  ml  hi  .  ai 

uses  =  1.  0  (4  for  trailing  edge),  and  the  location 
n 

for  may  be  left  blank  for  these  strips]  . 

If  alphas  do  not  vary  with  strips: 


1 


(12b) 


(13) 


Ui,  a^,  .  .  .  ,  aj^sZ 

If  alphas  vary  with  strips: 


(a) 

^1 1  ^2’  ■  ■ 

•  ’  ^ISZ 

(b) 

a^,  a^i  •  . 

•  •  °-isz 

(c)  cij,  * 

•  •  ’  ®ISZ 

V/b  w  series 
r 

(a)  (V/b^w)^, 

(V/b^w)^, 

numbe  r 


(b)  (V/b^a)j,  (V/b^w)^, 

numbe  r 


for  first  Mach  number 
for  second  Mach  number 


for  MSZ  Mach  number 


.  .  .  ,  first  Mach 

.  .  .  ,  (V/b  w)tc'7  second  Mach 
r  j  o 


(c)  (V/b^co)j,  (V/b^w)^. 
number 


(V/b^w)jg2  for  MSZ  Mach 


B.  Input  Data  Description 

(1),  (2)  Heading  Card  1  and  Heading  Card  2  may  contain  any  characters 
desired  in  Columns  2  through  72.  These  cards  are  convenient  for 
identifying  the  vehicle,  surface,  date,  engineer,  etc.  Both  cards 
may  be  blank  but  must  be  included  in  the  data  deck. 

(3)  Control  card:'  FORMAT  (1814) 

(a)  NTHRY  =  0,  piston  theory  is  used  to  compute  and 


25 


I 


NTHRY  4  0,  Van  Dyke's  theory  is  used  to  compute  and 
^2  A  =  0,  then  with  either  theory  and 

^2  are  the  same) 

(b)  NTHICK  =  0,  when  thickness  integrals  are  computed 
NTHICK  i  0,  when  thickness  integrals  are  given  (in  this 

case  they  are  constant  for  the  surface) 

(c)  NALPHA  =  1,  the  alphas  are  constant  (do  not  vary  with 

each  strip) 

NALPHA  =  ISZ,  the  alphas  vary  with  each  strip 

(d)  NTAUS  =  1,  the  r,  t,  ,  and  are  constant  for  all  strips 

h  t 

NTAUS  =  ISZ,  the  r,  Xj^,  and  x^.  vary  with  each  strip 

(e)  NZETAS  =1,  ^  and  L  are  constant  for  all  strips 

m  h 

NZETAS  =  ISZ,  and  vary  with  each  strip 

(4)  Control  card:  FORMAT  (1814) 

(a)  ISZ  =  number  of  strips,  <  25 

(b)  MSZ  =  number  of  Mach  numbers,  <15 

(c)  NO  PUNJ  =  0,  or  blank,  when  punched  card  output  is 

desired 

NO  PUNJ  0,  no  punched  output  is  desired 

(d)  JSZj  =  number  of  (V/b^a))'s  for  first  Mach  number,  <  20 

JSZ_  =  number  of  (V/b  («))'s  for  second  Mach  number, 

M  IT 

<20 

JSZ-,^_  =  number  of  (V/b  w)'8  for  last  Mach  number,  <  20 

MS  Z  r  — 

(5)  Single  parameters:  FORMAT  (6E12.8) 

(a)  sec  A,  secant  of  leading  edge  sweep  angle 

(b)  b  ,  reference  semichord 

r 
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(c)  s,  wing  semispan 

(d)  S,  wing  area 

(e)  c,  mean  aerodynamic  chord 

(6)  /^y.  series:  FORMAT  (6E12.8) 

Ayj^  .  .  .  ^yj32’  widths 

(7)  b.  series:  FORMAT  (6E12.8) 

’  ^ISZ’  semichords 

(8)  c^.  series:  FORMAT  (6E12.  8) 


'al 


.  c  TPr»i  control  surface  chords;  in  the  absence  of  a 
aioZ 


control  surface,  c^.  may  be  zero  or  blank,  but  a  sufficient 
number  of  cards  must  be  included 


(9) 


d^  series:  FORMAT  (6E12.8) 

di  .  .  .  distance  between  forward  and  aft  control  points 


(10)  Mach  number  series:  FORMAT  (6E12.8) 

Mach,  .  .  .  Mach. in  any  order  desired,  but  the  number 
1  MoZ  ' 

listed  must  agree  with  MSZ 
(11a)  Thickness  integrals  given;  FORMAT  (6E12.8) 

,  .  ,  .  ,  I^,  the  complete  airfoil  thickness 


(a) 

(b) 


h’  ' 

integrals 


Jj,  J2. 


.  ,  Jg,  the  control  surface  thickness 


integrals,  use  only  when  c  .  /  0 

clX 

^1’  ^h2’ 


,  .  .  .  ,  ^^i^^’^Bionless  chordwise 


coordinate  (Xj^/c )  for  the  control  surface  hinge  line 

(11b)  Thickness  integrals  are  computed:  FORMAT  (6E12.8) 

(a)  T^,  and  airfoil  thickness  ratios  (t/c)  at  point  of  max¬ 

imum  thickness,  hinge  line,  and  trailing  edge,  respectively 
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(b)  ^  .  and  dimensionless  chordwise  coordinates  for 

mi  hi 

point  of  maximum  thickness  and  hinge  line 

(12a)  Alphas  do  not  vary  with  strips  (alpha  is  a^,  the  initial  angle  of 
attack).  FORMAT  (6E12.8) 


number 


'^MSZ 


(degrees)  are  tabulated  in  order  for  each  Mach 


(12b)  Alphas  vary  with  strips;  FORMAT  (6E12.8) 


a^,  . 
number . 
(card). 


’  ‘  ’  °''ISZ  ®-^®  tabulated  for  each  Mach 

The  series  for  each  Mach  number  starts  on  a  new  line 


(13)  V/b^cj  series,  reference  reduced  velocity:  FORMAT  (6E12,8) 

There  is  a  reduced  velocity  series  for  each  Mach  number; 
each  series  starts  on  a  new  line  (card),  and  the  number 
of  V/b^oj's  must  agree  with  the  JSZ  for  the  respective 
Mach  number. 


C.  Example  Keypunch  Forms 

Example  keypunch  forms  are  given  on  the  following  pages.  Columns 
73  through  80  are  reserved  for  data  deck  identification.  This  space  may  be 
used  in  any  fashion;  however,  it  is  suggested  that  the  last  three  columns  be 
used  for  sequencing.  Only  the  cards  with  sequencing  in  Columns  73  through 
80  are  to  be  used  in  the  sample  data  deck;  the  lines  (cards)  with  Columns 
73  through  80  blank  are  for  clarification  of  input. 
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SECTION  IV 
PROGRAM  OUTPUT 


A.  Printed  Output 

1.  All  input  data 

2.  Thickness  integrals  (I's  and  J's) 

3.  Each  group  of  aerodynamics  influence  coefficients  (comprising 

a  complete  aerodynamic  matrix),  associated  Mach  number, 

and  V/b  co 
r 

4.  Sequencing  numbers  (Columns  73  through  80)  of  the  first  and 
last  punched  cards  (output)  for  each  group  (one  V/b^to)  of 
influence  coefficients 

5.  Example  problem  printed  output  is  shown  on  the  following  pages 
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B.  Punched  Output 

1.  A  deck  of  punched  cards  (output)  from  this  program  is  suitable 
as  an  input  deck  to  other  programs  requiring  the  use  of  AICs. 

2.  All  punched  output  is  sequenced  in  order  on  Columns  73  through 
80  starting  with  HMl  10000.  The  data  is  punched  in  the  following 
order: 

a.  Card  1  contains  and  M^:  FORMAT  (6E 12.  8) 

b.  .  Card  2  contains  the  size  (number  of  control  points)  of  the 

AIC  matrix  and  the  number  of  strips:.  FORMAT  (1814) 

c.  The  AIC  matrix  punched  in  column  binary  form  and  its  TRA 
card  make  up  the  remainder  of  the  punched  output  for 

(V/b^co)j 

3.  The  order  of  Statement  2  above  is  repeated  for  all  reduced 
velocities  and  associated  Mach  numbers  per  input  deck. 

4.  Each  AIC  matrix  is  punched  by  columns.  Column  1  starts  in 
Origin  1  and  Column  2  in  Location  (1  +  matrix  size). 

5.  The  oscillatory  AIC  matrix  is  punched  in  the  order  --  Column  1 
(real),  Column  1  (imaginary),  Column  2  (real).  Column  2 
(imaginary),  .  ,  .  ,  Column  N  (real).  Column  N  (imaginary). 

In  the  steady  case  all  columns  are  real  and  are  punched  in  order. 
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SECTION  V 


PROCESSING  INFORMATION 


Ope  ration 

STANDARD  FORTRAN  MONITOR  system 
Estimated  Machine  Time 
T  =  time  in  minutes 

ISZ  =  number  of  strips 

JSZM  =  total  number  of  reduced  velocities 

MSZ  =  number  of  Mach  numbers 

n  =  number  of  sets  (decks)  of  input  data 

T  =  1.  0  +  .  02  [(ISZ  •  MSZ  •  JSZM)^  +  (ISZ  •  MSZ  •  JSZM)^ 

+  •  •  •  +  (ISZ  •  MSZ  •  JSZM)^] 

Machine  Components  Used 

Core  storage,  about  5300 

Standard  FORTRAN  input  tape  (NT APE  2) 

Standard  FORTRAN  output  print  tape  (NT APE  3) 

Standard  FORTRAN  output  punch  tape  (NT APE  7) 


SECTION  VI 


PROGRAM  NOTES 


A.  Subroutines  Used 

RDLN,  reads  and  prints  title  cards 

AEROP4,  punch  AIC  matrix 

BINPU,  column  binary  punch 

All  other  subroutines  are  on  library  tapes 

B.  Generalized  Tapes 

Input,  print, and  punch  tapes  in  this  coding  are  defined  as  Units  2,  3, 
and  12,  respectively;  however,  these  maybe  altered  by  placing  the  desired 
units  on  symbolic  cards  HM110060,  HM110061,  and  HM110062. 
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FLOW  DIAGRAM 
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SECTION  VIII 


SYMBOLIC  LISTING 


Some  of  the  symbols  used  in  the  program  are  defined  as  follows: 


FORTRAN  Symbols 

NTHRY 

NTHICK 

N  ALPHA. 

NTAUS 

NZETAS”^ 

NO  PUNJ 
.  ISZ 
MSZ 

J  SIZE  (M) 

JSZ 

SEC  LAM 

BR 

S 

CAP  S 
C  BAR 


Definition 

Option- -theory  used  for  C^,  C^ 

Option- -thickne  s s  integrals  given  or 
computed 

Option- -q's  constant  or  vary 

Option--T's  constant  or  vary 

Option--|'s  constant  or  vary 

Option- -punching  or  no  punching 

Number  of  strips 

Number  of  Mach  numbers 

Number  of  reduced  velocities  for 
Mach  number 

Number  of  reduced  velocities  for  a 
Mach  number 

sec  A 

b 

r 

s 

S 

c 


Please,  no  remarks  about  our  Greek! 
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SYMBOLIC  LISTING  (continued) 


FORTRAN  Symbols 

Definition 

C  BAR  I 

C  BAR  2 

^2 

RAD  DEG 

w/180.  0  (program  constant) 

DELTA  Y(I) 

Ay  for  strip  i 

B  (I) 

b  for  strip  i 

CA(I) 

c  for  strip  i 
a 

D  (1) 

d  for  strip  i 

EMACH  (M) 

m'th  Mach  number 

EKR  (J,  M) 

l/k^  =  (V/b^w)  for  reduced  velocity  j, 
for  m'th  Mach  number 

El  (N) 

I  series  (thickness  integrals) 

EJ  (N) 

J  series  (thickness  integrals) 

AI  (I.  N) 

I  series  for  strip  i 

AJ  (I.  N) 

J  series  for  strip  i 

ZETA  H  (I) 

for  strip  i 

ZETA  M  (I) 

£  for  strip  i 

’m  ^ 

TAU  (I) 

T  for  strip  i 

TAU  H  (I) 

T,  for  strip  i 
n 

TAU  T  (I) 

for  strip  i 

ALPHA  (I,  M) 

a  for  strip  i,  for  m'th  Mach  numbe 

EK  (I,  N) 

K  series  for  strip  i 
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SYMBOLIC  LISTING  (continued) 


FORTRAN  Symbols 

Definition 

CQJNST  (I) 

2 

4(b/b^)  Ay/s  for  strip  i 

A  (I.  N.  K) 

Premultiplying  matrix  in  oscillatory 

coefficients  matrix  equation 

G  (N,  K) 

Real,  oscillatory  leading  edge 

coefficient  matrix 

Gr(N.  K) 

Imaginary  matrix 

H  (N.  K) 

Postmultiplying  matrix  in  oscillatory 
coefficients  matrix  equation 

Q  (N.  K) 

Working  array 

QI  (N,  K) 

Working  array 

P  (N,  K) 

AIC  matrix,  complex 

The  symbolic  listing  of  the  program  is  shown  on  the  following  pages. 
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HM11C002 


o 

o 

o 


GO  O' 

o  o  o 
o  o  o 
o  o  o 


Z  z  X 
z  z  z 


i|  fO 


o  o  o 
I  o  o  o 


z  z  z 
z  z  z 


u> 

tu 


l/)ll 

(/>  m 

<NJ 

^  H  •  I 

^  A 

<  o  ^ 

i  H-  fNi  — 

H  Z 

o 

i/>  in 
<N  r-* 

--  II 

o  z 

*/M! 


^  m 

sO  <M 

*  ^ 
in  Q  < 

ni  iu  z 
<«»  X  a 

X  ^  ^ 

z  < 


^  •  m 
^  ^ 

«  * 
tn  <-•  o 

N  •  <NI 

||«  n  z 

||<  «-  X 

•  Ui 

z  z 
^  a 

ae  V) 
o  z 


z  a 

11.^2 


n  < 

UJ  H 

•  UJ 
Nl 

o  «> 

•  K 
O 

•  D 
< 

a.  K 

•  * 
a.  z 

X  3 

•  < 

•-•  H 

o  - 

•>  3 
O  < 


fsl  O 
^  «» 
<n  < 
u 

n  • 

*  A 

X  • 
U  > 
< 

z  < 

h» 

UJ  ^ 

•  UJ 

<  o 
z  • 

a  z 

-J  UJ 

^  • 

ac  < 

z  « 

UJ  n 

< 
z  « 
o  < 
z  • 
z  « 
a  UJ 
u  *1 


11  o  I 

I  ni  rg 
|o  o  o 
o  o  o 


z  z  z 
X  z  z 


o 

u 

UJ 

o 

z 

Ui 

3 

3 

u. 

z 


< 

z 
> 

o  >• 
o  oc 
ac  o  I 

Ui  Ui 

<  z 
z  ^ 

in 

<n  z  , 
O 
♦  ^  I 

X  m  I 
.  fn  ^ 
GO  a» 

• 

nj  >► 

I  O  A  I 
llu.*  Z  Z 
pO  'O 


<  < 
z  z 

«  cc 
o  o 

rg  m 


In  'O 
rg  rg  fM 
o  o  o 
o  o  o 


z  z  z 

Z  X  X 


llz  •X 
_i  A  ^ 

II 3  • 
a  A 
UJ  ^ 
i||A  UJ  A 

.  r-4  t 

o 

Q  II 
Ui  UJ 

o  <  ^ 

3  O 
O  A 

III  ^  ” 

Hloc  < 

3  n 
jx  I 
A  K  fn 

ll  ^  z 
<  » 

“OX 
UJ  A 

[  hr  n  n 

I  X 

11^  n 
^  z 


X  v. 
00  ^ 

A  «r  •• 
UJ  A 
AO* 
Z  A 


Z 

o 

<  H 
Z 

3  A 
Z  <  A 
m  ^  I 
h*  O  “T 


!|M  X 

^  UJ  n 

•  o 

X  X  z 
^  ^ 
'T  ^  X 


•H  rg  <n 

fn  ^  fO 

o  o  o 
o  o  o 


XZZl 


UJ 

^  Xl 
fn 


n 


-•  z 


X  X 


z 

A 


X| 
X 
X  X  I 


X 

A 


X 

a.  r“* 

A  & 
K  •-• 
l/>  A 
Z 

A  00 
«««  ««4  Z 

n 

» 

♦-*  X 
^  X  I 

A  in  I 

* 

»-*  o 
nj  X  o  ll 

UJ  ^  z 


A  A 

o  o 


<0 
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1E16.8  }  )  HM110037 

FORMAT  (IH  IIX,  29HALPHA  ZERO  SERIES  (DEGREES)  =  6(1F8.2>  3X  )  HM11003e 


HM110039 


W-  ITN  nO 
\-t  ^  >t 

o  o  o 
i  o  o  o 


z:  s:  s: 

XXX 


o  ^  Cvj 
in  u> 
o  o  o 
o  o  o 


X  X  £ 
XXX 
00 

•»■«*  • 
•>  <i 
OO 


M|^0  K  00 

Mm  lA  in 
|o  o  o 
Mo  o  o 


z  X  X  <: 
X  X  X|: 


nj  m 't  |j 
UO  'C  sO  II 
o  o  o  I- 
!o  o  o  I- 


X  X  X  ' 
IIX  X  X; 


loo  O'  o 
l-o  vO  r- 
o  o  o 
o  o  o 


i'x  X  X 
I'X  I  X 


o  o  o 
o  o  o 


XXX 

I  X  X 


*0 

</>  I 

,  Ol  l( 


«  X 

X  m 
fM 

i/)  *-4  » 

-J  O' 
<  ^ 

a:  <si 

o 

LU  “n 
H  X  ^ 

z  ^  ^ 


CO  X  ^  I 
to  rsj  J 
UJ  “5  1 
Z  X 
X  ^  I 
O  ^  f 

X  “)  X  ! 

X  Ai 

'I* 

o 

tu  •  — 
X  m 
O  — 
a 

X  a  X 

o  ^  < 

U  X 
X  »-  • 
to  X 
<N  X  A* 
m  ^ 

m 

X  •* 
m  X 
4*  vn  >»* 

|o  O  -J 
I  X  X  X 


--  fli  X 

w  o  a  %  ; 

(/)  <  Ml*  M  : 

I  M  X  X  ,  jg 

i  X  m  X  III  I 

X  ^  O'  iW  M 

l  ml  1 


X  V  X  HI  jjui  ^ 

i£  <  I  B|j|  jj|  ®  0 

to  '*  X  1  II  I  tu  L 

I  <  ^  o  m  «  I  AJ  f 

^  <  K  M  I  ^  t 

>-  X  -J  X  : 

o  ^  X  • 

<  ..M  jW 

uj  X  o  nniiiiiiitif  I  * 

A.  to  H  l\  CD  « 

to  rH  o  K  Kb  I  • 

X  X  111  i\-o ' 

mh  X  -H  Mu  ' 

I  o  '  I  ' 

I  •  m  •  • 

lx  X  X  • 

fA  ^  \\m\\\  X  : 

m  X  m  I  jan  m  < 

o  o  I 

X  X  ■  X 


<  <  1  Hi  <  <  < 

X  z  1 1  z  X  X 

a  OC  ■H  OC  OL,  CC 

o  o  'll  KPl p  o  o  o 

ii  I 

■o  BHii  o  Al 

'i  ^  i  H  i  AI  AJ  rg 


I  eg  m 
•-«  o  M 
]\  H  O 
to  Uj 
I  at  o  o 
a  OC 

<  <  o 
t-  o  < 

Z  Z  OC 


to  A4 
a.  to 

►  <  I 

t*  - 

z  ^ 

o  •  » 

Q.  to  •-< 

O  • 

Z 

OC 
*  CO 

eg  >- 

</>•»< 
X  X  H 
<  ^ 
M  ^  ilj 

eg  o  o 
to  ui  *• 
•«i  to 


AJ  AJ  eg 
at  at  at 

Q.  0.  Ou 
<  <  < 


at  at  UJ 

CL  (L  ^ 

<  <  < 


K  H  A 
=>=>=? 
Q.  a  CL 
Z  Z  z 


o  o  o 
<  <  < 
UJ  UJ  at 
oc  QC  a£  I 


3 

a 

^  o 

o  o  I 

<  o 

UJ  O 

OC  “> 
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CO  202  I=1»MS2  HM110077 


rsi  ff)  >t 
QO  ao  00 
o  o  o 
o  o  o 


GO  9^  O 
00  CD 

o  o  o 
o  o  o 


u>  -o 

O'  O'  ^ 

o  o  o 

O  O  Q 


o  f-i  ru 
o  o  o 


'O  oo 
o  o  o 


Z  X  z 

X  Z  X 


ill 


z  z  z 

z  z  z 


z  z  z 

X  X  z 


O  -J  I 
O  iu 
o  a  Ni 


i/)  H 

Z  M 

«k  * 

•  of 
(O  A 

««  o  o 


M  uj 

««  3 

<  z 

o  ^ 

^  H* 
N  Z 
iL  <  O 

M  z  U 


3  A 

a  ^ 

Z  fx  A 
N 

o  o 

<  ^ 

Ul  o  3 
oe  o  o 


f>J  o 

A  N 
a  (M 
<  » 
o 
Z  N 
N  Nj 
UJ  (/) 

Q.  ^ 


I  H 

K  A 
3  < 
a  K 
Z  A 
«  NJ  fV 
Z  fVJ 

c  ^ 

< 

Ui  U.  O 

ec  ^  q\ 


rsi 

r>»  NI 
1/)  #■* 


<-•  P-4  o 

II  < 
I  «  ^ 
A  11 

3  ^ 

<  in 
^  (M  «« 

z  m  3 
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226  00  228  I^ltlSZ  HM110112 

IF  ICAIIJJ  227.227,228  HM110113 

227  ZETAHm-1.  HH110114 


At  1,3,31-Bin/CAm 
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EKi  I,1J  =  U.0/E  MACHIM}  )  •(CBAR1+2.0*CBAR2*E  MACHJ  M)  •AJ  (  I,  1 )  HM110226 

♦3.0*CBAR3*EMS*IAni,4)+ALPHA|I,M)*ALPHAnfM)  )  )  HM110227 

EK<  I,2)  =  (  1.0/E  MACHIM)  )«lCBARl-t^4.0>CBAR2»E  MACHtM)  HM110228 


I 


O'  o 

III 

ir\  o 

H 

■'in 

111  rg  fO  ^ 

eg  ff) 

H  {jjlHI  III! 

HI  fn 

>»•  «f  -I 

Ii  ^  ^  mi 

1  liiil  in  in  in 

<N  eg 

11  IllfU  III 

rg  nj  rg 

rg  rg  rg 

Hi  rg  rg  rg  liil 

yw  rg  rg  rg 

o  o 

D  lllo  III 

o  o  o 

o  o  o'p  « 

o  O  o  i 

g  o  o  o 

n  IlliU  lly 

pH  ^  ^ 

pH  pH  pH  tJ^JN 

pH  pH  pH  ^ 

llllli  pH  pH  pH 

II  llinil  III 

pH  pH  pH 

pH  pH  pH  «4  4H 

^  pH  pH  m 

1  llll  pH  pH  pH 

z  z 

1  H  B 

z  z  z 

■  ii )  II  ■  ij 

Z  Z  Z  IWH 

Z  Z  Z  Q 

i  X  z  z 

I  X 

z  X  z 

""  lii  li:i 

XXX  ■  1 

XXX  1 

I  X  X  X 

i  ll 

lii 

fO 

II  III!  Ill' 

lii 

yiy 

|ll|!llllllj 

yiy 

Hill !!! 

ll|ll|lll|1| 

Z  w 

nllH 

mil!!**! 

Hiiiii|m 

•  M 

illll  mil  !!.!! 

iiiiiiiiiiii 

< 

1 

miiii|im 

1 

« 

llilillllllllll 

llll|l||l|y 

mill|l!!y 

X  z 

iimiiim 

a  w 

I  I 

ni  z 

■pH  HHnnl 

1 

-i  x 

1 

oc  • 

Ill  III 

^  PH  P«S 

1 

<  u 

m||ll| 

—  X  imiuiiiii 

1  i 

*  < 

ml|ll| 

A  w 

X  «»  1 

^  z 

|l||||!ll 

O  <. 

III  llll  ii  il 

PH  M  1  1  II  1 

z 

•  X 

4  w  «p  m{|j{[f||i 

1  j 

•  UJ 

yiiiiii 

1  *  ^ 

<  i[[ii[iiiiii 

llllllll] 

r 

f-4  <1 

»llliii 

rg  ^ 

!l!  llll  III! 

ui  <  X  iiiiiiimii 

ljilll 

1 

^  rg 

mIII!I| 

♦  < 

!!!  Illll  ,1!' 

rg  «  a.  Illllll[|||| 

1 

1 

<  oc 

yiiim 

Ill  I'll!  Ill 

1  _ 1  [HiHiiiiii 

miiij 

1 

X  ^ 

mjiii 

PH 

o  z  <  ill 

ijiiii 

a. 

ii  mil  iiii 

fjHilill 

illl  rg 

-J  o 

yiiiiH 

ilHl 

bO 

<  « 

mijiii 

««  M 

pH 

♦  o 

m||ll| 

X  — 

m  •  '  II'II 

t 

mi 

<  “J 

1 

p 

in  sO 

lllilll 

H-  < 

•  ♦ 

UJ  — 

ill  'ill'  iif' 

illll 

Z 

pH 

ym 

fsl  « 

III  |:::|  ||||| 

a  ^4  X  BO 

llllli  HP 

—  OC 

i 

nil 

ii^ 

X 

•-I  < 

O  Z  PH 

llll  K  III! 

Him 

lj 

<  (D 

y||i|M 

•  UJ  PH 

IIIIIIIIIIII  < 

«  O 

Mil  II'II  7|| 

mi 

Ii 

llliillll  z 

o  ^ 

yjllly 

— 

1 

UJ 

•  « 

•  OC  PH 

1.  ; 

liiimmij 

rg 

1  1 

^  PH 

III;  illll  III 

^  •» 

n  1 

oe  CO 

\f\  II 

rn  ^ 

♦ 

<  O  X 

1  ^  in! 

rg  -H 

00  Z 

llli|y 

A  « 

II  liii:  ill;; 

Z  T  BO 

il 

M 

z  — 

llllllj 

u  o  < 

^^^K^lllll  llll 

p  ffi 

11  *D  H 

UJ  X 

1 

PH  •  H- 

II  Illll  IS 

1 1 

mi 

r*  r> 

•  o 

^lllljj 

•  ^  UJ 

UJ  ii 

rg  UJ 

m  <t 

yjllllj 

PH  ♦  fSj 

»  & 

oc  z 

yjllllj 

^  1 

Z  Z  m  ■ 

i^BaSS: ' 

<0  < 

< 

y|||j|j 

PH  pH  O 

■  t  II  Mil 

msa  y 

hp  >- 

CD  UJ 

y|||jij 

Z  t>  « 

in  '111'  i'"j 

EKSl  S 

rg  z 

o  « 

yjllllj 

M  PH  «4 

llliillll  Illll 

•  o 

ymilj 

i  w  w 

llli  '!!;:  IS 

O  UI  z  I 

HI  ii 

UJ 

o  « 

Qlllj 

o  ^  • 

•  Nl  UJ  II 

1  ^  ^ 

1  ^  a. 

• 

yjllilj 

<  <  PH 

iSi  Illll  1 1 

[k^w^n  iiiiiiiiiii 

III  1  ^  fflrj 

< 

yjllllj 

z  «  z 

lllllll 

tllllllll  '*?  W  Im 

PH  h* 

♦  'V 

y||lllj 

Ul  PH  • 

injillll 

iiiiiiiii  «k  mil 

Z 

«<»  o 

millj 

Z  ^ 

mi!J||j]j 

llliillll  PH  pH  8»  11m 

ik  K 

rg  8 

liijiij 

«  HP 

llll  liil-  I'-l 

.!'*mitii 

IllillZ  N  3  i 

->  ^ 

•  'T 

imilj 

X  < 

u  IS  M 

II  lllllll! 

llliillll*^  *9  a.  BU 

PH  0. 

miiij 

—  ox 

'•lllllll 

llliillll  UJ  ^  IP  E 

OC^  00 

—  M 

imiij 

M  <  a 

mil Hll 

lU^X^  lll'!l|ll 

1  ^  ^  9ni  f 

X  3  N- 

•-•  PH 

PH  X  -J 

UJ  8  PM  g 

III  1  PH  (D  n  all  f 

UJ  O  cv 

<  fO 

imiij 

^  UJ  < 

<0  '•J  |»|  ■  I 

llliillll^  o  nil 

•  • 

umij 

»  «  « 

II  1  1 

llliillll'^  <n  UJ  lUI 

PH  UJ 

•h« 

Ml 

PH 

9gi  r 

iiiiiiiii  K  ^  Km  ! 

H-  o 

mJIIII 

1)11111  1  IhJ  PH  M 

PH  M 

Mj 

X 

Z  B 

ii|v)  O  Z  1 

U.  z  O 

UJ 

mjP 

UJ 

•u  i 

*7  O  X  ■ 

1  z  u 

pH 

pH  rg 

llll 

^EEl  II  H 

Ii 

n 

yjllllj 

pH 

HHH 

1  in  ^ 

1  rr  r* 

1 

^  ill 

[  II 

■ 

rg  rg 

«  9*  O 

I  m  in  nC 

fsi  f>i  rsj 

o  o  o 


ac  z  z 

XXX 


GD 
O 

►  t  v» 

I  o  o  oc 

MHO 

I  ^  ^  QC 

•»  •»  S 

1  *v  M 

V9  O  O) 

o 

CO 

tM 
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282  G( 1,3)>-EK< 1»A}*BB  HM110264 

GI2,3)*-EK( I,5)*B8  HM110265 

GI3>l)-0.0  HM110266 


304  WRITE  OUTPUT  TAPE  NTAPE3.  21.  ( ( P( K. L > I ) . L^l . IN) , K= I , NU)  HM110340 

30^  CONTINUE  HM110341 

IF  I  NOPUNJ  I  308.306.308  HH110342 


I 


f-4  O  tA 
^  O  A| 

^  >o 


^  <0  O' 

«o  rg 

(A  00  CO 

(S4  ^  ^ 

<A  ^ 
O  Q. 
UJ 


o  I 

o 

UJ 

z 

UJ  u 
^  UJ 

<  o 


o  ^  M 

^  ^O  tA  I 

^  4*  ^  fl 

4*  4*  4- 

o  o  o 

<0  AJ  OD  i 
A  lA  ^  1- 

|m  m  m 

|AJ  AJ  CM 


O  KJ  T 
O  A  2 
O  £  D 
a. 
o 
z 


iH  rg  ^ 

A  A  'O 

^  o  r* 

^  fA 

o  o  o 


O  <A 
|0  4*  AJ 
1  At  ^  Al 
[A  ^  A 


O  h-  O' 
^6  4-  f- 
fA  00  fA 


I  “>  •<  «  I 
iU  O  i 


z  u 

«  o 

o 


A  A  i 
xA  4)  A  I 
4  4  4 
4  4  4 
O  O  O 

h«  lA  ^ 
A  A  4 
(A  fA  <A 
AJ  Al  AJ 


H  Z  A 
>-•  Q 
X 
< 
u 


z 

< 

A  X 
fA  O 
4  O 
4  X 
O  O. 

(A  UJ 
<A  O 
rA  X 
AJ  ^ 
O 
A 
A 

^  OC 

o 


AJ  Al  O 
O  A*  f'- 
AJ  O  h- 
4  4  rA 
O  O  O 


^  *.4  (V 

A  ^  ^  H 


^  <A  ^ 

A  eo  f*- 
fA  00  eo 


M  X  UJ 

Ul  O  AJ 

< 

Z  A 

Ul  -7 


A  O 
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case,  [fJ  =  (l/2)pV^S/c)  jh]  .  The 

piston  theory  is  limited  to  high  Mach  number  (or 
high  reduced  frequency),  but  Van  Dyke's  quasi¬ 
steady  correction  extends  the  validity  to  some 
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